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POWER TRANSFER TO ION CYCLOTRON WAVES 
IN A TWO ION SPECIES PLASMA 
by Donald R. Sigman and John J. Reinmann 
Lewis Research Center 


SUMMARY 


Calculations were made on the efficiency of using a Stix coil to transfer power to 
waves near the ion cyclotron frequencies of a cold, collisionless plasma consisting of two 
ion species. The mass ratio m, /m of the ion species was 2, and the frequency region 
near and between the cyclotron frequencies (f2- a and of each ion species was inves- 
tigated. For wave frequencies oo less than 52^, only the usual ion cyclotron resonance 
absorption peak is obtained. For J2 ia > co > two absorption peaks are obtained. 

The additional peak corresponds to power absorbed by fast wave modes associated with 
the heavier species. The wave set up nearer the cyclotron frequency of the lighter ion 
S2 ia is a left-polarized wave, and power is absorbed primarily by modes whose wave- 
lengths, which are parallel to the static magnetic field, are near the wavelength of the 
Stix coil. The second absorption peak occurs at lower frequencies where the power is 
absorbed by the right-polarized fast wave modes associated with the heavier species. 
There is a sharp cutoff on the low frequency side of this second absorption peak at 


co = 


°ia n ib 


Xfl ia + Sb 
XS2 ib + n ia 


1/2 


where x is the density ratio of the two species, n^/n • When x is made small the 
latter peak disappears. 


INTRODUCTION 

Previously, extensive calculations were made (ref. 1) on the efficiency of trans- 
ferring power from a Stix coil to waves in a cold, collision-free hydrogen plasma at 
frequencies near the ion cyclotron frequency. In these calculations, the plasma was 
assumed to contain only electrons and a single ion species, and the effects of electron 
density, inhomogeneities, Stix coil wavelength, etc., on power transfer were investigated. 



However, plasmas used for the study of ion cyclotron waves often contain more than a 4 
single ion species, such as H + , Hg, D + , and D*. Extending the calculations of refer- 
ence 1 to include the effects of a second ion species therefore has practical as well as 
theoretical interest. Since the mass ratios D + : H + , and D^D 4 " are all 2 and 

since it would be unusual to have a plasma consisting primarily of and H* ions, a 
mass ratio of 2 was chosen for these calculations. In this report calculations are pre- 
sented to show how power transfer varies with magnetic field in the region near and 
between the cyclotron frequencies of the two ion species as the relative concentration of 
the two species is varied. 


THEORY 


The basic equations to be solved are the same regardless of the number of species 
present in the plasma. These equations are given in references 1 and 2 where Maxwell's 
equations and the particle equations of motion (in polar or rectangular coordinates) are 
combined and written for a cold, collisionless plasma as 


V X V X E = 


tj- K • E 
cr — 


( 1 ) 


where it has been assumed that all field quantities are of the form 

F(r,t) = 


and 



( 2 ) 


(3) 


The terms S, P, and D in the dielectric tensor K are given in chapter 1 of reference 2. 
They depend on the wave frequency oj , the plasma frequencies 


n 


2 

k 


w 2 

m k 


(k = e,a, b) 


( 4 ) 


of the various species, and the cyclotron frequencies 
2 



(k = e,a,b) 


(5) 



where k is the wave number in the direction of the static magnetic field. The coef- 
ficients aj and are determiaed by the boundary conditions. 

The term P is determined primarily by the electrons, so that the terms S and D 
are the ones most affected by the presence of a second ion species. The variations of v j 
and i >2 with frequency are determined mostly by the properties of S and D. The 
terms u 1 and v 2 are given by 



2 2 2 

Figure 1 contains plots of and against k for several frequency regions. In 
some of these plots, there are two branches for which k^ is positive. A branch which 
intercepts the k^ axis at points marked R has right polarization, while those having 
intercepts at points marked L are left polarized. The species a and b are labeled 
so that m./m =2. In this report, only those plasmas are considered whose electron 

u -14 o 

density is so low (usually below 10 cm” ) that the right polarized wave is not able to 

propagate in a finite radius plasma column ■tfhen only a single ion species is present. 

2 2 2 

This right polarized wave has k < (u /c )R. 

First, consider the region in the frequency spectrum (fig. 1) to the right of w = 
where is the cyclotron frequency of the heavier mass ion species. The calculations 
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show that here the presence of species a only changes the power transfer characteristics 
from a Stix coil to the ion cyclotron wave quantitatively and not to a large enough degree 
to elaborate on in this report. 

Because the density is low, there can be no wave propagation at all for frequencies 
to the left of = Q. , where is the cyclotron frequency of the lighter species. 

Between the frequencies f^ a and (fig. 1) there are two special points marked 
D = 0 and S =0. The terms S and D can be zero in the frequency range 
S2ia> u) > ft ib only when more than one ion species is present; it is the way S and D 
vary with frequency and/or magnetic field that leads to some interesting power transfer 
results. The conditions for S = 0 are 



^ia n ib! 


X fi., 
a lb 


X £2. 

a ia 


lVia\ 

+ ¥ib/ 


(8a) 


and for D = 0 


2 2 Y 0 2 . 
w = W D ” X b^ia + 


X ft?, 
a lb 


(8b) 


where X & and X^ are the relative concentrations of the species (X & + X^ = 1). Waves 
which can be set up by the Stix coil in the region to the left of D = 0 (co > co^) are left 
polarized and those to the right are right polarized. From equation (1), it can be shown 

O 

that Eg/E r = iD/(n - S), and as has been pointed out by another investigator (private 
communication with M. A. Rothman of Princeton University) Eg = 0 at D = 0 so that 
the Stix coil should be ineffective in coupling power to plasma waves at this point. 


RESULTS AND DISCUSSION 

Calculations of power transfer were made from a Stix coil to waves in the plasma 
for o» > ft^ with the relative concentrations of the two ion species the principal variable. 
For these calculations, the Stix coil wavelength was 45 centimeters, the coil radius 
10 centimeters, the plasma radius 5 centimeters, and electron as well as the total ion 

13 

density was 10 per cubic centimeter. The results of these calculations are shown in 

figure 2 where relative power transfer P* (defined in ref. 1) is plotted against 

ft = o>/ft. for several X. /X n ratios. The term P* is the total power transferred to 

la. Da. \ O s|tO 

the plasma wave from the Stix coil divided by 2na j 
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where j is the amplitude of the current density in the current sheet (i.e. , Stix coil) 
which extends from -a < z < a. 

As the relative concentration of the heavier species is increased, the following 
features can be seen in figure 2: 

(1) The amplitude of the resonance peak closest to u = decreases. 

(2) The resonance peak moves closer to w = and becomes narrower. 

la. 

(3) A second peak appears to the low frequency side of the first resonance peak. 

(4) The amplitude of the second peak increases and the peak becomes narrower. 

(5) The resonance peak of the heavier species near co = moves to the right and 
broadens, as expected (ref. 1). 

Calculations show that for the resonance peak nearest f2 ia , power is absorbed 
primarily by modes whose wavelengths, which are parallel to the static magnetic field, 
are near the wavelength of the Stix coil and that the wave is left polarized. This 
resonance peak is the normal ion cyclotron resonance peak associated with the lighter 
species. To show how the heavier species affects the position of this resonance peak, 
the arrows in figure 2 indicate the position of the peak if the heavier species were absent 
and the ion and electron density were equal to the ion density of the lighter species in the 
mixture. For the absorption peak for which o> D > u> > co g , the wavelength of the modes 
absorbing the power is considerably longer than the Stix coil wavelength, and the waves 
are right polarized. 

In addition, it should be noticed that there exists a sharp cutoff on the low frequency 

side of the peak just to the left of w = u> s , and, for the conditions of the model, no power 

is transferred to the plasma between this cutoff point and u> = Also, the amplitude 

of the resonance peak near u> = is not severely reduced until the ion density of the 

ID 19 

heavier species becomes less than about 10 per cubic centimeter (ref. 1). 

Thus, the presence of the lighter species does not severly affect power transfer 
near the cyclotron frequency of the heavier species, but the presence of the heavier 
species can severely affect power transfer near the cyclotron frequency of the lighter 
species. 

2 2 

An examination of the curves for v against k (fig. 1) reveals that the modes 
absorbing power in the peak nearest cu = o>g are the fast wave modes (i. e. , waves with 
lower k values) associated with the heavier species. When the lighter species is 
absent, these modes absorb power at different frequencies, and there are distinct 
absorption peaks for each of the modes. When a lighter species is introduced, all the 
fast wave modes associated with the heavier species are squeezed together, and all 
absorb power near w = u>g. In addition, there are fast wave modes associated with the 
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lighter species which absorb power for u > O ia - The absorption peaks for these modes f 
are not shown in figure 2. 

To examine the effect of increasing the electron density, a power transfer calcula- 

14 

tion was made in the region SL a > u> > for n g = 10 per cubic centimeter, X a = 0. 4, 
X, = 0. 6, and all other conditions the same as those for figure 2. This power transfer 

b i4 

curve is given in figure 3, where it can be seen that, at n g = 10 per cubic centimeter, 
the ion cyclotron peak has virtually disappeared, while the peak near w = u>g is much 
larger. This occurs because the fast wave modes associated with the heavier species 
now have wavelengths close to the wavelength of the Stix coil when co is near cog. For 
n g < 10*2 per cubic centimeter, there is no peak near co = Wg, since the fast wave modes 
do not absorb power efficiently at low densities. 


CONCLUSION 

When a second ion species is included in the calculations of power transfer to 
plasma waves in the vicinity of either ion cyclotron frequency, additional features are 
found in the absorption spectrum. As the relative concentration of the heavier species 
is increased, a peak, which corresponds to power absorption by right polarized waves 
associated with the heavier ion species, appears near u>g (wave frequency at S = 0). 

In addition, the normal left polarized cyclotron peaks are changed in magnitude and 
position. For high densities (n g > 10 14 ), this peak, corresponding to right polarized 
waves, dominates the spectrum. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 31, 1967, 

129-01-05-09-22. 
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APPENDIX - SYMBOLS 


a 

one -half of Stix coil length 

r 

radial cylindrical coordinate 

B 

]r 2 (P + S) u> 2 /PS + S 2 - D 2 \ 

S 

defined in ref. 1 



S c 2 \ S / 

t 

time 

B 0 

static magnetic field parallel to X 

ratio of heavier to lighter ion 


z -axis 


species concentration 

C 

P 

k 4 - 2 — k 2 S + — (S 2 - D 2 ) 

X 

relative ion species concentra- 


S 

2 4 


tion 



c c J 



c 

speed of light 

z 

axial cylindrical coordinate 

D 

defined in ref. 1 

Z k 

charge number for particles of 





type k 

E 

electric field vector (time 




varying) 

e 

azimuthal cylindrical coordinate 

e 

electronic charge 

u l’ v 2 

defined by eq. (7) 

F(r,t) 

electric or magnetic field 

n 

plasma frequency 


quantity 

0 

cyclotron frequency (or fre- 

f(r) 

radiaLvariation of F(r,t) 


quency ratio) 

.* 

J 

CO 

current density of wave exciting 

wave frequency 


current sheet 

COg, u> D 

u> when S or D equals zero 

k 

wave number 

Subscripts: 

K 

dielectric tensor 

a 

lighter ion species 

L 

defined in ref. 1 

b 

heavier ion species 

m 

species mass 

e 

electron 

n 

density 

i 

ion 

P 

defined in ref. 1 

k 

particle species 

P* 

relative power transfer to 

r 

radial coordinate 


plasma 

z 

axial coordinate 

R 

defined in ref. 1 

e 

azimuthal coordinate 
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(c) Concentration ratio, X a :X b , 0.6: 0.4. 

Figure 2 . - Relative power transfer as function of fre- 
quency ratio (Q = Stixcoil wavelength, 

45 centimeters; coif radius, 10 centimeters; plasma 
radius, 5 centimeters; electron number density, 
10 13 per cubic centimeter. 
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Relative power transfer to plasma, P , ohms 


3. OxlO 4 



Figure 3. - Relative power transfer as function of frequency ratio 
for electron number density of Ur 4 per cubic centimeter. Stix 
coil wavelength, 45 centimeters; coil radius, 10 centimeters, ■ 
plasma radius, 5centimeters; concentration ratio, X a ;Xjj, 

0.4; 0.6. 
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